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CN , Abstract 

^ I In previous studies, we have examined a resonant excitation of disk oscillations 

! in deformed disks. In these studies, however, mathematical treatment around the 

I resonant points was not rigorous. In this paper the inadequate point is corrected, 

""I. with no essential changes in the final results. For this excitation process to work, 
disks must be general relativistic. That is, the non-monotonic radial distribution of 

O ■ epicyclic frequency in relativistic disks is essential for the presence of the resonance 

^ I and for trapping of oscillations. In this paper, the growth rate of resonant oscillations 

li^l is expressed in a form more suitable for numerical calculations. 

^ i Key words: accretion, accretion disks — black holes — high-frequency quasi- 

^ ! periodic oscillations — neutron stars — relativity — stabihty — X-rays; stars 

^ ! 
(N ■ 

oo ■ 1. Introduction 

o ; 
oo . 

O ■ In previous papers (Kato 2004, 2008a), we proposed a resonant excitation mechanism 

^ ■ of disk oscillations in deformed disks. The purpose is to propose an excitation mechanism of 
^ ' quasi-periodic oscillations (QPOs) observed in neutron-star and black-hole X-ray binaries (e.g., 
d I Kato and Fukue 2007, Kato 2008b). In this resonant excitation model of disk oscillations, a 
deformation of disks from an axially-symmetric equilibrium state is essential. The deformation 
to be considered is a warp or an eccentric deformation of disks in the equatorial plane. 

An outline of the model is as follows. A non-linear coupling between a disk oscilla- 
tion (hereafter we call it the original oscillation) and a deformed part of the disks (warp or 
eccentric deformation) brings about some forced disk oscillations (we call them intermediate 
oscillations). The intermediate oscillations make a resonant coupling with the unperturbed 
disk at particular radii of the disk. After this resonant coupling, the intermediate oscillations 
feedback to the original oscillation by a nonlinear coupling with the deformed part of the disk. 
Since this nonlinear feedback process involves a resonance, the original oscillation is amplified 
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or dampened. Kato (2004, 2008a) examined this nonlinear feedback processes, and derived an 
excitation criterion and growth (damping) rate of the resonant oscillations. It is of importance 
to note that in Keplerian disks this resonant excitation process works only when the disks 
are general relativistic. That is, a non-monotonic radial distribution of epicyclic frequency is 
necessary for appearance of resonance and for trapping of oscillations. 

In these studies the intermediate disk oscillations were assumed to be local near to the 
resonant radius in the sense that their radial wavelength is shorter than the characteristic 
radial length of disks. That is, the radial derivative, d/dr, operated to the wave quantities 
was taken to be ik. Further, k was assumed to be constant. This treatment oi k = const., 
however, was inadequate around the resonant point. If this inadequate treatment is properly 
corrected, we find that unlike in the previous papers, the resonant point is not the place where 
a local dispersion relation of the intermediate oscillations is satisfied with a constant /c, but the 
radii of Lindblad resonance for intermediate oscillations (when resonance occurs by horizontal 
motions) or the radii where the intermediate oscillations are trapped in the vertical direction 
(when resonance occurs by vertical motions). 

In this paper, we modify the analyses in our previous papers so that the above-mentioned 
inadequate treatment is corrected. We apply the mathematical techniques used by Meyer- 
Vernet and Sicardy (1987) in their study of resonant disk-satellite interaction. The results 
show that the stability criterion derived in the previous papers are unchanged. The expression 
for growth rate obtained in the previous papers is found to be still applicable, but in this paper 
the expression is changed in a form more suitable to numerical calculations. 

We first summarize in section 2 the basic equations and relations necessary to the study 
of the present resonant excitation problem, although they are given in the previous papers. In 
section 3, the stability criterion and growth rate are derived to the case of pressure-less disks as 
it is instructive, and to the case of disks with pressure in section 4. The final section is devoted 
to a brief discussion on a meaning of the instability criterion. 

2. Brief Summary of Basic Equations 

We briefly summarize here basic equations and relations [see Kato (2008a) for details] 
that are necessary in the subsequent sections. 

2.1. Nonlinear Hydrodynamical Equations 

In the present problem, the general relativity is essential. For simplicity, however, for- 
mulation in this paper is done in the framework of a pseudo- Newtonian using the gravitational 
potential introduced by Paczyhski and Wiita (1980), assuming that the central object has no 
rotation. We adopt a Lagrangian formulation by Lynden-Bell and Ostriker (1967). 

The unperturbed disk is in a steady state with a steady flow Uq. By using a dis- 
placement vector, ^, a weakly nonlinear hydrodynamical equation describing adiabatic, non- 
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self-gravitating perturbations is written as, after generalizing the linear equation derived by 
Lynden-Bell and Ostriker (1967), 

PO0 + 2p,{u, ■ V)^ + m = poC(t 0, (1) 

where L{$,) is a linear Hermitian operator with respect to $, (Lynden-Bell and Ostriker 1967) 
and is 

L{^) = Poiuo ■ V)(uo ■ V)^ + po(^ ■ V)(V^o) + V [(1 - rOpodiv^ 

-poV(div^) - V[(| ■ V)po] + (I • V)(Vpo), (2) 

and po{r) and po{r) are the density and pressure in the unperturbed state, and Fi is the 
barotropic index specifying the linear part of the relation between Lagrangian variations 6p and 
6p, i.e., (5p/po)iincar = Ti ((5p/po)imear- Siucc, the sclf-gravity of the disk gas has been neglected, 
the gravitational potential, V'o(^); is a given function and there is no Eulerian perturbation of 
it. At the stage of equation (1), there is no restriction on the form of Uq, but in the main part 
of this paper, no is taken to be a pure cylindrical rotation. 

The right-hand side of wave equation (1) represents the weakly nonlinear terms. 
Derivation of a detailed expression for C is boring, because of the presence of many terms. 
However, after lengthy manipulations we can summarize them into a relatively simple form in 
the case of Fi = 1: 

PoC(^,0 = -^M.e.,^(V^o) - T^fpol^VeA (3) 



2 dvidrj dvi \ drj 

using Cartesican cordinates. When Fi ^ 1, some terms should be added on the right-hand side, 
which are shown in equation (82) of Kato (2008a). An important characteristics of C is that 
we have commutative relations (Kato 2008a) for an arbitrary set of ^3- 

/ Po^M^2.^3)dV = J Po^iC{t,,^,)dV = J p,i^C{i,X2)dV. (4) 

As shown later, the presence of these commutative relations is the reason why we have a 
simple stability criterion. We suppose that this commutative relation is a general property of 
conservative systems. 

2.2. Nonlinear Coupling and Growth Rate 

We now assume that the disks are deformed from an axisymmetric steady state by some 
external or internal cause. The deformation is assumed, for simplicity, to be time-independent, 
and has the azimuthal wavenumber mw- The Lagrangian displacement associated with the 
deformation, ^^(r,t), is denoted by^ 

i (r,t) =exp(-imw</^)C {r,z). (5) 

^ Here and hereafter, ^ [sec equation (6)], and [see equation (7)] are expressed in complex forms. 
When their nonhnear coupHngs are calculated, we must be careful since the real displacement associated 
with deformations, for example, is not ^ itself, but the real part of ^ given here. 
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As the deformation, a warp will be the most probable, but it is not only one candidate of 
possible deformations. A plane-symmetric one-armed spiral deformation is one of another 
possible candidates that can excite disk oscillations. In both cases we have mw = 1. 

Our purpose here is to examine how the behavior of disk oscillations are affected by 
the disk deformation. As shown below, some of oscillation modes are resonantly excited on 
deformed disks through nonlinear coupling with disk deformation. The nonlinear coupling 
processes are schematically shown in figure 1 of Kato (2004). 

The displacement vector associated with a disk oscillation, ^, is assumed to have fre- 
quency uj and azimuthal wavenumber m. Separating the time and azimuthal dependences from 
^(r,t), we introduce ^ as 

^(r,t) = exp[i(ci;t — m(y9)]^(r, 2;). (6) 

The first step of the nonlinear interaction between the disk oscillation characterized 
by {(jj^m) and the deformation characterized by (0,1) introduces two kinds of intermediate 
oscillations with azimuthal wavenumber being m+1 and m — 1. Let us denote the displacement 
vector associated with these intermediate oscillations by 

^l\r,t)=eMi{^t-^^)]Q{r.z). (7) 

where m = m-|-l orm — 1. Here, represents the intermediate oscillations resulting from the 
coupling between ^ and ^ , while does those resulting from the coupling between ^ and 
^ , where the asterisk represents the complex conjugate. That is, and ^„ are described, 
respectively, by 

-uj^P^C + 2«^Po(^xo ■ V)^+* + L{^1') = \ [poC{lf)+poC{f,i)l (8) 

-cVor + 2^c.po(t^o ■ V)r + L{t) = \[poC{lf*) + poC{f*M. (9) 
Next, the second stage of the nonlinear coupling is considered, which is a feedback 

^ int 

process returning to the original oscillation, ^, by the intermediate oscillations ^_|_ interacting 

- W - int 

with ^ . In the case where the coupling occurs through , the feedback is described by 

-u'p,i + 2iujp,{u,.V)k + L{l) = ^[poCCC\t*)+PoCit\C)]- (10) 

^ int 

On the other hand, in the case where the feedback occurs through ^_ , the equation corre- 
sponding to equation (10) is 

-uj'poi + 2zu;po{uo-V)k + m = l[poC{Ct) + PoC{t,C)]. (H) 

An important point to be noted here is that as a result of this feedback process, the 
original disk oscillation is amplified or dampened, i.e., the frequency, u, can be no longer real, 
since in the feedback processes a resonance is involved. How much is the imaginary part of a;? 
This can be examined from equations (10) and (11) using the fact that the operators ipo{uQ ■ V) 
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and L are Hermitian (Lynden-Bell and Ostriker 1967). After some calculations the imaginary 
part of frequency, tUj, can be formally expressed as (e.g., Kato 2008a) 

= (12) 

where ± denotes the cases of the coupling through and respectively, and 

= f 3/ \p.r\c(C\i"") + c(^"":C)]'iv, (13) 
= y^/ ^''»«'ic(«-'.«") + c(f:c')]dv, (14) 



E = \^„l Po€ ^ - i(m ■ VMdV, (15) 



and 

E = , 
2 

where c^o is the frequency of the original oscillation before the mode couplings. The above 
expressions for W± and E have physical meanings that W± is the rate at which work is done 
on the original oscillations by the nonlinear resonant processes, and E is the wave energy of 
the original oscillations [see equation (93) by Kato (2001)]. 

It is important to change the above expressions for W± by using the commutative relation 

(4) as 

W- = ^^J lpot[C{l\t) + CC^',k*)]dV. (17) 

2.3. Oscillations in Isothermal Disks 

Hereafter, we restrict our attention to oscillations in geometrically thin disks. The 
steady unperturbed disks are axially-symmetric and have no motion except for rotation. Here, 
cylindrical coordinates (r, ip, z) are employed, in which the 2;-axis is perpendicular to the disk 
plane and the origin of the coordinates is at the disk center. The unperturbed flow is then 
described as no = (0,rr2(r),0), where f2(r) is the angular velocity of disk rotation. We further 
assume, for simplicity, that the disk is isothermal in the vertical direction. In the vertically 
isothermal disks, the unperturbed density in disks, po{r,z), is stratified as (e.g., Kato et al. 
1998) 

r -^^ 

Po(r,^) =Poo('^)exp 



2m{r)V ^^^^ 
where poo is the density on the equatorial plane, and H is the half-thickness of the disk and is 
related to the vertical epicyclic frequency, fi^, by 

nlH' = ^-^ = cl{r\ (19) 
Po 

where Cs is the isothermal acoustic speed. The vertical epicyclic frequency, fi^, is equal to 
the angular velocity of the Keplerian rotation, fix? in the case of the central object being 
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non-rotating, and practically equal to il, since the disk is assumed to be geometrically thin. 
Hereafter, however, we use Q± without using fix or Q, so that we can trace back the effects of 
Q± on the final results. 

Concerning oscillations superposed on such disks, we consider oscillations whose radial 
wavelength is moderately short so that the radial variations of physical quantities in the un- 
perturbed state [including the radial variation of H{r)] are neglected compared with the radial 
variation of wave quantities. Furthermore, we assume that the oscillations also occur isother- 
mally. Then, we can neglect the term of V[(l — ri)podiv^] in the operator L in equation 
(2). 

If the above simplifications and approximations are adopted, the r- and z- dependences 
of ^{r,z) are approximately separated. That is, the 2;-dependence of $,{r,z) is obtained by 
solving an eigen-value problem, and the eigen-f unctions are found to be an orthogonal set of 
Hermite polynomials (Okazaki et al. 1987). Restricting our attention to one of such oscillations, 
we write ^(r, z) as 

i{r,z)=^r,n{r)nniz/H) (20) 

Ur,z)=^^,^ir)nniz/H), (21) 

Ur,z)=iAr)nn-i{z/H). (22) 

where Tin is the Hermite polynomial of argument z/H, and n{= 0,1,2,...) characterizes the 
number of node(s) of oscillations in the vertical direction. It is noted that the number of 
node(s) of in the vertical direction is smaller than those of C,r and by one, as shown 
in equation (22). However, the subscript n (not n — 1) is attached to as ^z,n in order to 
emphasize that ^r,n, C^.n, and ^z. set of solutions. 

Under these approximations, we express the r-, ip-, and z- components of the homoge- 
neous parts of wave equation (1) as 

r r^^ 1 V w 

- (iu - mn? + - An"^ - c"^— ^rn-i2n(uj-mn)^^n + ^lH^^ = 0, (23) 

dr^i ' dr 

-{uj-mny^^ + i2n{LO-mn)Cr = 0, (24) 



[-{LU-mny + nnl%^r,-nnlH^^ = 0. (25) 

Nonlinear coupling between ^ and introduce intermediate oscillations The 
azimuthal wavenumber of the intermediate oscillations are m + 1 or m — 1. Their z-dependence 
are characterized by n + 1 or n — 1 when riw = 1 since TiiHn = Ti-n+i + nHn-i and by n when 
nw = 0. To consider these various coupling cases separately, we write ^_|_ in the forms of 

i^{r,z)=C%d^)n^(^/H), (26) 

i$X{r,z)=^l%^^{r)nn{z/H), (27) 
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(28) 



where the subscript n represents n + 1 or n — 1 or n. 

Nonhnear couphng terms are also separated into terms proportional to exp[— i(m± 
and Tinir/z). That is, in the case of coupling through we write the coupling terms as 

\po[C{i,e) + C{eMr = PoY.K+,n{r)^M^{^t-rhvWn{z/H) + ■■■ (29) 



Po[C(^,r ) + C(r ,^)]^ = Po E^^,+,-(0exp[2(^t - m^)]'Hn{z/H) + 



(30) 



Po[C(^,r ) + CiCMz = PoJ2^z,+Ar)eMK^t " mif)]nn-i{z / H) + ■■-, (31) 



where H — denotes terms orthogonal both to Hn and Hn-i, and the subscript + is added to A's 
in order to emphasize that they are related to the 93-dependence of exp[—i{m + l)ip]. In a similar 
way, in the case of coupling through the intermediate oscillations of (l/2)po[C'(^,^^*) + 
C(^^*,^)] can be expressed in forms similar to equations (29) - (31), introducing A^.-.n, A^,-,n, 
A 

The equations describing intermediate oscillations are then written as 

^2 



dr 



A, 



fint 



dr 



A 



(32) 
(33) 

(34) 



where m represents m + 1 or m — 1 , and n does n + 1 or n — 1 when nw = 1 and n = n when 
nw = 0. From equations (32) and (33) we can eliminate C^X,n 



JZcmt 



-A 



(35) 



dr ' ' u — rhQ 

Hereafter, we use equations (32) - (35) as the basic equations describing the intermediate 
oscillations. 

Finally, we should notice that in the case of isothermal disks described above, W.^. and 
W- given by equations (16) and (17) are written in the forms of 



and 



UJQ 

2 



^ j p^,{x)[{2-Kf''n\rH]- 
[ Poo{x)[{27Tf'fi\rH]- 



(^int A* I fint A 



"int A* I i*int 4* I j^* 



1 

n 



dx. 



dx. 



(36) 
(37) 
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where [( comes from the part of volume integration in the Lp- and 2;-directions. 

Furthermore, the wave energy, given by equation (15) can be expressed as 

E = ^^u;l{r^Hpoo)cEn, (38) 



where En is a dimensionless quantity given by 

|2 



2 



E„ 



TU — ^ "''^-+ "--1 !^ — , 39) 

and the subscript c denotes the values at resonant radius, which is defined by Ji = or J2 = 
(see below). It is noted that in the case of = 0, (n — 1)! is zero. 

3. Resonant Excitation of Oscillations in Pressure-less Disks 

Before examining general cases of disks with pressure, we study here the limiting case 
of pressure-less disks, i.e., Cg = and H = 0. Equations (35) and (34) show that ^^^^ (and ^™*) 
and respond resonantly, respectively, at the radii where —(u — mf2)^ + = holds and at 
the radii where —{u — rfiQ)^ + nQ\_ = holds. We define here Ji{r) and J2('") by 

j^l^r) = -{00 - rhny + k'^ and J2{r) = -{u - rhn)^ + nQl. (40) 

We call the resonance at Ji = horizontal resonances (Lindblad resonances) and the resonance 
at J2 = vertical resonances. We consider these two resonances separately. 

As a typical example in infinitesimally thin disks, we consider here the case where all 
motions including those associated with disk deformations are in the equatorial plane, i.e., 
n = 0, riw = 0, and h = 0.^ At the radii of Ji = 0, and resonantly respond to A's as [see 
equations (35) and (33)] 

. = i f - ^-^T,A,,^,^) , (41) 



Ji V ' ' CO — mf2 



/•int 



1 

Jl 



2Q ^ -{uj-mQ)^ + K,'^-4ri 



2 



(42) 



uj — rhD. ' ' {uj — rhiiy 

The quantity ^f*" can be taken to be zero, since it does not respond resonantly and has no 
contribution to the final results. 

To examine the resonant interaction at the radii of Ji = 0, we introduce a small imaginary 
part of u, i.e., Ui, a.s u = 00,. + iui, where the imaginary part, Ui, is tentatively assumed to 
be negative so that causality is satisfied. Later, the results obtained by c^i < is extended 
analytically to the whole region of Ui as usually done in stability analyses. Hereafter, the real 
part, Ur, is written by u or Uq, without confusion. 



^ This expression for En is different from that given in the previous paper (Kato 2008a), since the previous 

one is not suitable to represent the sign of wave energy. 
^ In the foUowing formulations, however, we retain n, nw, and n in general without specifying to particular 

values. 
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Near to the resonance, Ji is expanded as 



Ji = J[ 



.2{uj — mQ)c 



-UJ; 



where the subscript c represents the values at the resonant radii, and J[ is 



m(LU 



■mQ)Q 



dlnQ 



,dlnK 



(43) 



(44) 



dr dr 

Let us now consider an integration of /(r)/ Ji, / being an arbitrary smooth real function, along 
the radial direction including the resonant region. The imaginary part of the integration comes 
from the path near to the pole. Thus, we have 



fir) 
Mr) 



■dr 



— tt/c/ J'l when {uj — mQ) / J( > 
^/c/ J'l when {uj — rfiQ) / J( < 0. 



This can be summarized as 
fir. 



rn 

J Jii 



-dr 



fc 

vr— ^sign(u; 



rfiQ)^ 



(45) 



(46) 



where sign(ti; — mQ)c is the sign of ci; — friQ at the resonant radii. 

Let us now estimate the rate of work done by the resonance on the original oscillations 
by using equation (46). In the case of coupling through the work is given by equations 
(36) and (37). Since .^"^ are given by equations (41) and (42), respectively, we have 

from equation (36) and (37) 



W± = -^7r(27r 



N3/2 



^,irHp^ 



n\- 



00 J c 



1^(1 



sign(ci; — mVL)^ 



A. 



2n 



r,±,n 



— l- 



UJ - 



-A 



(47) 



where we have used /^exp(2;^ /2H^)nl{z/H)dz = {2TT)^/^h\H, and subscript H is attached to 
Ui in order to emphasize that this is a case of horizontal resonances. Hence, the growth rate, 
—Ui [see equation (12)], can be expressed as 



vrn! 



-sign(tt; — mf2)( 



Ar 



2n 



mil) 



A 



(48) 



2rl\J[\uJoEr, 

It is noted that the above expression for growth rate is valid even when n ^ 0. Equations (47) 
- (48) formally hold even in the case with pressure, as shown in the next section. 

The vertical resonance might be unrealistic in pressure-less disks, but we briefly sum- 
marize the formal results here, since the expression for growth rate is applicable even in the 
case of disks with pressure, as shown in the next section. We now expand J2 around J2 = as 



J2 — J'o 



2{u-mn)c 

[r - rc) - I uJi 

J'? 



where 
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miuj ■ 



mQ)Q 



dlnQ 

dr 



dlnQj 
dr 



(49) 



(50) 



and Tc is the resonant radius of J2 = 0. In the above two cases of the horizontal and vertical 
resonances, the resonant radii are different each other. However, the same notation, r^., is used 



here and hereafter without confusion. 

After this preparation, we calculate the rate of work done on oscillations by the same 
procedures as those in the case of horizontal resonance. The results show that the growth rate 
is given by 



Tvih — iy. 

-^i,Y,±,h = -TTTTiri E^sign(t^ - mQ), 

Resonances in Disks with Pressure 



A 



2 



(51) 



Equations describing the intermediate oscillations [equations (32) - (34) or equations 
(35) and (34)] are now solved without assuming Cg = and H = 0. Since we are interested in 
behaviors of ^ in a region close to the resonant point of Ji = or J2 = 0, all the coefficients of 
^ in equations (32) - (34) are assumed to be constant except for Ji or J2. That is, in the case 
of the horizontal resonance (Lindblad resonance), J2 = const., and Ji is taken as [see equations 
(43) and (44)] 

Ji = J[x — 2i{uj — 'rhQ)cUJi, (52) 

where x is the radial distance from the resonant radius, Tc, i.e., x = r — re- In the case of the 
vertical resonance, on the other hand, Ji = const., and J2 is taken as [see equations (49) and 
(50)] 

J2 = — 2i{uj — mQ)cUJi, (53) 

where x = r — r^. The resonant radii in the horizontal and vertical resonances are different each 
other, but the same notation, r^, has been adopted here without confusion. 

In disks with pressure, different from the case of pressure-less disks, the resonant region 
is not narrow. That is, the resonant region is widened by the effects of pressure (see below), 
and thus the coupling terms, A^^ A^, and A^, can be no longer regarded to be spatially constant 
in general. Considering these situations, we introduce Fourier transforms of ^ and A as 

|"*(A;) = / exp{-ikx)^ {x)dx, (54) 



' —00 

00 



A{k) = / exp{—ikx) A{x)dx . (55) 



After these preparations, to solve the set of equations (32) - (34), we apply a method of Fourier 
transform adopted by Meyer- Vernet and Sicardy (1987) in their study of resonant disk-satellite 
interaction. Here and hereafter, such subscripts as ± and h to be attached to ^ , A, and 
A are sometimes omitted in order to avoid unnecessary complications, unless they are explicitly 
necessary. 

4.I. Horizontal Resonances (Lindblad Resonances) 

In the present case, the Fourier transform of equations (35) and (34) give, respectively. 
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drh 



and 



(56) 



(57) 



where e = —2{ijj — ihVL)uJi. 

Eliminating ^"*(A;) from equation (56) and (57), we have an inhomogeneous differential 
equation oiS,^^^{k) with respect to k: 



dk ' J[ 
where A{k) is defined by 



+ -fjC\k) 



2Vt - 
uo — mil 



.Cgk'^ (i^ — friVl)'' 
J'l J2 



kn'^.H - 

i^^A,. 

J'?. 



Let us first consider the case of e/J[ > 0. Equation (58) is then solved to lead 



e-*(A;) = A+)(fc)exp 



J'l 3J[J2 



where 



1 



A{k')exp 



Ji 3J1J2 



dk'. 



(58) 



(59) 



(60) 



(61) 



In determining the integration range of equation (61), we have used a boundary condition 
allowing .^"*(A;) for k —00. This boundary condition is required, since C'^'(/c) is the 
Fourier transform of ^™*(x). 

By using the inverse Fourier transform: 



^f{x) = — / CJ'\k)exp{tkx)dk 



271 



we can express C'^*(x) in an integration form as 



27r 

where ap is defined by 



e;'\x) = — Z^(+)(A;)exp 



—k + i[ kx- ^alk^ 
■Ji 



1 

—I 

3 



dk. 



Cg(a; — mVl)' 
J( J2 



Since i^^^{x) and ^™*(x) are related to by [see equations (33) and (34)] 



2fi 



A, 



[to — miiy 



S2 



1 

J2 



int 



dr 



A., 



(62) 



(63) 



(64) 



(65) 



int 



the work integral, VFh,+ , in the case of the couphng through is written as, using equation 
(36), 



Wn,+= y53 J dx J^^dk^poo[{27Tf^^n\rH]D^+){k)exp 



—k + i[ kx - ^apP 
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X 



A* 



^* 



[X) 



—ikVl?, HA* 
J2 



[X] 



(66) 



where real parts on the right-hand side of equation (66) have been omitted. 

In the case of e/ J[ < 0, the expression for .^™*(A;) is shghtly changed from equation (60) 

as 



er(fc) = /^(_)(A;)exp 



-k 



where D(^^){k) is given by 



yl(A;')exp 



:k' 



3J{J2 



k' 



k 



>3 



dk' 



(67) 



(68) 



.J'l 3J{J2 

where A{k) is defined by equation (59). The integration range in equation (68) is changed from 
equation (61) so that the boundary condition = is satisfied at k = 00. By performing the 
same procedures as those in the case of e/J{ > 0, we can derive an expression for 1Vh,+. An 
expression for Wu,+ in this case is the same as equation (66), except that is now changed 



to D 



(-)■ 



I.e., 



2 



dx 



dk—p„[^2„f'^h\rH]D^.^(k)e^p 

D ZTT 

2Q 



J'l 



k + i{ kx - -a%k^ 



(x) + —iknlHAl 
■J2 



[X] 



(69) 



int 



Expressions of the work integral in the case of the coupling through ^_ ^ can be derived 

^ int 

by similar procedures as those in the case of the coupling through ^. In this case, the 
subscript + attached to A's in equations (66) and (69) is changed to — and m = m — 1. 
4.I.I. The case of constant coupling terms 

First, we consider the case where the coupling terms, Ar, A^, and Az, can be regarded 
to be spatially constant in the resonant region. (This is really realized in a particular case, as 
mentioned in the next section.) An another purpose of examining the case of constant coupling 
terms is to have a rough image on the width of resonant region in disks with pressure. 

If Ar, Aip, and Az are spatially constant, their Fourier transforms give 



Ar{k) = 27iArS{k), A^ = 27rA^S{k), A^ik) = 27rAzS{k). 
In this case, D(j^\{k) given by equation (61) is reduced simply to 



2ti 

D^+){k) = -i—H^{k) 



UJ 



2Q ^ 



(70) 



(71) 



where H+{k) is the unit-step function, i.e., H^{k) = for /c < 0, while H{k) = 1 for k> 0. Since 
D(^^){k) does not involve x, the integration with respect to x in equation (66) leads to 27r6{k). 
Hence, performing the integration with respect to k in equation (66), we have 



-^7r(27r^3/2 



(ri^'poo)c 

nl 



J'l 



2n . 2 



(72) 
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in the case of e/ J( > 0, where m = m + 1 and n = n — 1, or n, or n + 1, depending on the type 
of couphngs and the form of deformation. 

In the case of e/J[ < 0, is given by equation (67). By performing the same 

procedures as those in the case of e/J[ > 0, we can derive an expression for PFh,+- The final 
expression for Wi{,+ in this case has the opposite sign from equation (72). 

Inequahty e/ J[>0 imphes that {u — mfl) / J[>0, since —Ui should be taken tentatively to 
be negative as mentioned before, while inequality e/ Ji < means {u — rfiQ) / J( > 0. Considering 
this difference, we can summarize 1Vh,+ in the above two cases of e/J[> and e/ J( < as 



^0 /o \3/2 
TT 27r ' 

2 ^ ' 



nl — ■ — - — sign(ct; — m\L)^ 



A 



r,+,ri ' 



UJ 



2Q ^ 



(73) 



where m = m + 1 . 

In the case of coupling through we have, after the similar processes as the above. 



-— 7r(27r^ 



3/2-1 V HOO c . I 

' n\- — , sign(u; 



1^(1 



mf2)( 



UJ 



2VL ^ 



(74) 



where m = m — 1 . 

Finally, from VFH,±,fi given above and the wave energy, given by equation (15), we 
find that the growth rate, equation (38), is written as 



vrn! 



-sign(ti; — mf2)( 



20. 



-A 



2r-^\Ji\ujQhn uo — mil 

which is identical to equation (48), and also equal to the growth rate derived in the previous 
paper (Kato 2008a), as will be mentioned at the end of this section. 
4.1.2. The case of radially changing coupling terms 

The arguments in subsection 4.1.1 show that in the case where Ar, A^, and Az are 
radially constant, the radial structure of resonant region is determined by [see equation (66) 
and notice that D(^^-^{k) has the step-function] 



3? / dk exp 



jk + iikx- ^a^k^ 



(76) 



where the term with e can be practically neglected. The a;-dependence of this integral has been 
examined by Meyer- Vernet and Sicardy (1987) [see equation (41) and Fig, 5 in their paper]. 
They show that the half width of the resonant region is ~ ap, which is on the order of (rH'^Y^^, 
i.e., r(H/ry^^, in the present case. 

In general, Ar, A^, and A^ radially oscillate with the wavelength determined by the 
product of ^ and ^ . Hence, in some cases the scale-length of the variation of A's will be shorter 
than r(if/r)^/'^, although it depends on where the resonance occurs (see the next section). To 
examine general cases, numerical calculations are needed. However, in order to have a rough 
image of the effects of spatial variation of Ar, A^, and A^ on the work integral M^H,±,ri; we 
consider here a limiting case where Ar, A^, and A^ oscillate periodically in the radial direction. 
That is, we assume 
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Ar{x) = A^^''exp{ikox), A^{x) = A^^'^exp{ikox), A^{x) = A^^'^exp{ikox), (77) 

where A^^^'s and ko are the amphtudes and wavenumber of the spatial variations. Then, the 
Fourier transform defined by equation (55) show that 

Ar{k) = 27ri(°)5(A; - ko), A^{k) = 2nA^^'>6{k - ko), A,{k) = 27rif 5(A; - A;o).(78) 

Let us consider the case of e/J( > 0. In this case, i5(+)(A;) given by equation (61) is 
reduced to [cf., equation (71)] 



D(+)(/c) = -z^i7+(A;-A;o) 



exp 



e 

A' 



—ko+iYko 



.(79) 



Since Ar, A^, and A^ are given by equation (77), performance of integration with respect 
to X in equation (66) gives 27r6{k — k^). Next, we perform the integration with respect to k. 
Since D(+)(fc) is now given by equation (79), we have 

2VL 



2 ^ ' J[ 



i(o)_,MMi(o)' 
UJ — mil ^ Jo 



(80) 



Similar arguments can be made in the case of e/J( < 0. In this case, Wh,+ has the sign 
opposite to equation (80). Expressions for VFh,+ in two cases of e/J( > and e/J( < can be 
combined into a single form as before. Similarly, we can consider the case where the coupling 
occurs through ^_ . All of these results can be summarized as 



00 c 



\Ji\ 



sign(ct; — mQ)^ 



X 



A 



(0) 

r,±,n 



2^ 



— i- 



ijj — mQ 



koQlH 



-A 



(0) 

2,±,n 



(81) 



Hence, the growth rate is written as 

vrfi! 



2rl\J[\u^En 



sign(ci; — mVL)^ 



X 



4 (0) , 21] ^ (0) . kQn\H ^ (0) 



(82) 



uo — iJiQ ' J2 
^.2. Vertical Resonances 

In the case of the vertical resonance, from equations (35) and (34), we have, respectively. 



A 



and 



A.. 



(83) 



(84) 



2 dk 

Eliminating from the above two equations, we have an inhomogeneous differential equation 
with respect to as 



dl 



int 



ikn^l , H 



Ji 



Ar 



2Q 



A^]+A, 



(85) 



14 



where 



p,v 



JUi 



(86) 



In deriving equation (85), in order to avoid unnecessary complications, we have neglected /c^Cg 
compared with Ji, since the magnitude of kH contributing to the present formulations is smaller 
than unity, i.e., kH ^ 1. 

In the case of e/J2 > 0, equation (85) is solved as 



e(fc) = ^v,(+)(fc)exp 



e .1 3 g' 

^, rv 1~(y,ry \rrh 



J'. 



where the coefficient Dv,(+)(^) is given by 



1 r 

Dw,{+){k) =-i-j- J_ 



1 

2 

xexp 



ik'miH 
Ji 



jrk' + t-al^^k'^ 



Ar-i 

dk'. 



(87) 



2VL ^ 
uj — rhfl 



By using C"*(/c) given above, we can write and ^"^^ as [see equations (35) and (33)], 



Sr 



— Z- 



and 



tint 



Jl Jl 

2n A;fi2 H^.^^ 



2n - 

UJ — mil 



2VL 



2Q ^ 
UJ — rnVt 



UJ — mf2 Jl ^ Jiuj — mVt 
Substitution of these relations into the work integral (36) gives 

Wy,+= / dx r dk^poo[{27rf/'n\rH]Dy,(+^{k)exp 



—k + i\kx 

■J2 



X 



ikVt^ H / ^ 



Jl 



A* + i 



2^ 
UJ — rfiQ 



1 

-A* 

n 



(88) 



(89) 



(90) 



(91) 



In the case of e/ J2 < 0, Wy^+ has a similar form, but we should use, instead of Dv,{+)(^), 



Dy^(^^){k), which is given by 



J2 -^oo 



ik'nQ] H f - 

-^{Ar 

■Jl 



2VL - 
UJ — mil) 



A, 



xexp 



e 

1'. 



k' + i^alyk'^ 



dk'. 



(92) 



.1 

'2 3 

As in the subsection 4.1.1, we consider two cases where i) A^ is constant in the resonant 
region and ii) is spatially periodic as A^°)exp(i/cox). In the former case, summing the two 
cases of e/ J2 > and e/ J2, we have 



— -■n{2'Ky' {n-l)\ — — — sign(u;-mi2)c|v4^,+,n| • 

J < 



(93) 



and in the latter case we have 



^0 n3/2/~ t n, (''-f^Poo)c . 

■,n=-Y7r(27r) / (n-1)! — — — sign 



'uj — mQ)^ 
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X 



Jl 



2VL 



-A 



(0) 



A 



z,+,h 



(94) 



LU — m " 

In the coupling through the work integral, Wy^^^n, has the same expression as 
equations (93) and (94), except that the subscript + is now changed to — . 



After the above considerations, we summarize the growth rate, Ui 



lV,±,n: 



case of Az being constant as 

7T(n 



-sign(ci; — ml]); 



A 



In the case where A^ changes as A^ = A^^''exTp{ikox), we have 

7r(n-l)! 



in the former 



(95) 



X 



2r3|J^|cuoS„ 



signi^cj 



A. 



(0) 



2VL 



-A 



(0) 



A 



(0) 



Jl V ' ' oj — m ^' ' / ' ' 

Finally, we note here that the growth rates obtained in the previous papers are the 
same as those obtained in subsection 4.1.1, although the final expression in the previous papers 
are different from those in subsection 4.1.1. In the case of horizontal resonance, for example, 
Gh,± in Kato (2008a) is equal to rcJ(/2 in the present paper. Furthermore, if and are 
eliminated from equations (46) - (48) in Kato (2008a), we have 

J2 



2n ^ 

~M} '''' 



(97) 



D\ u- 

where D is J1J2 — clk'^{uj — rhiiy. This expression for shows that Cr defined by equation 
(52) in Kato (2008a) is written as 

2n 



Cr = ^2 A 



LU ■ 



(98) 



nQ']_)c into equation (70) 



Then, substituting Gh = fcJi/'2 and equation (98) with (J2)c = {—1^ 
in Kato (2008a), we obtain an expression for growth rate in the case of horizontal resonances, 
which is identical with equation (75) in subsection 4.1.1. 

Similar arguments easily show that the growth rate given by Kato (2008a) for vertical 
resonances, i.e., equation (71) in Kato (2008a), is identical with that in section 4.2, i.e., equation 
(94). 

5. Discussion 

In this paper we have derived a stability criterion on resonant excitation of disk oscilla- 
tions in a deformed disk. Deformation of disks that can lead to excitation of disk oscillations 
are a warp and an eccentric deformation of disk plane, i.e., rri^ = 1, although this is not dis- 
cussed in this paper (see Kato 2004, 2008a). In previous papers (Kato 2004, 2008a), the above 
resonant excitation problem was examined and a condition of excitation of disk oscillations 
was derived. From mathematical point of view, however, the treatment of disk oscillations 
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around the resonant point was not suitable in the previous papers. In this paper, we corrected 
the inadequate points. The results of the correction show that the stability condition and the 
growth rate derived in the previous papers are kept qualitatively unchanged. In this paper, 
however, the growth rate of oscillations is written in a form applicable more easily to numerical 
calculations of growth rate. 

In disks with pressure, the resonant region is widened by pressure effects. Based on 
arguments by Meyer- Vernet and Sicardy (1987), we showed in section 4 that the half-width 
of the resonant region is ~ r(/7/r)^/^ (or ~ H(r/Hy^^), which is shorter than r but longer 
than H. The coupling terms comes from nonlinear products of original oscillation, ^, and disk 
deformation, ^ . If the disk deformation is assumed to be global, the scale of spatial variation 
of the coupling terms is mainly governed by the spatial behavior of the original oscillation, 
^. Hence, in general, the width of resonant region is comparable with the wavelength of the 
original oscillations and may not be narrow enough to regard the coupling terms A's as constant 
in the region. Hence, in section 4, we examined two limiting cases where i) the coupling terms 
are constant in the resonant region and ii) they vary sinusoidally in the region. Analyses in 
section 4 show that in these two limiting cases there are no essential differences in growth rate, 
especially no difference in the criterion of stability. To examine what happens in intermediate 
cases of the above two limiting ones, however, numerical calculations will be needed. 

It is noted that the case where the coupling terms. A, are roughly constant in the 
resonant region is realized when resonance occurs near to the boundary between the propagation 
and evanescent regions of the the original oscillations. This is because near to the boundary, 
the radial wavelength of the oscillations is long. Applications of the present disk oscillation 
model to high-frequency QPOs by Kato and Fukue (2006) and Kato (2008a) are concerned in 
such cases. The another limiting case where the coupling terms vary in the radial direction 
will be realized when resonance occurs in the propagation region far apart from the boundary 
between the propagation and evanescent regions. 

There are two kinds of resonances, i.e., horizontal and vertical resonances. In both 
resonances, the condition of excitation of oscillations is that the wave energy of the original 
oscillation. En, and lu — friil at the resonant radius have the opposite signs, i.e., ' sign(cL' — 
mfl)c < 0. For the resonance to occur efficiently, the resonant radius must be within the 
region where both original and intermediate oscillations exists predominantly. Hence, it will be 
reasonable to suppose that the sign of {u — mQ)^ roughly represents the sign of the energy of 
the intermediate oscillation. [Notice that the sign of wave energy of the intermediate oscillation 
is determined by the sign of — rfiQ in the region where the wave exists predominantly, see 
equations (38) and (39)]. In other words, the condition of i5„ ■sign(u; — mf2)c < is roughly equal 
to the condition that the original and intermediate oscillations have energies of the opposite 
signs. This instability criterion thus can be understand as a result of energy exchange between 
the original and intermediate oscillations at the resonant point (Ferreira and Ogilvie 2008). 
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The resonant instability condition, i.e., En ■ sign(ci; — fhVL)c < 0, however, is not exactly 
equal to the condition of both oscillaions (original and intermediate one) having the opposite 
signs. This suggests that the energy exchange at the resonant point is not only between the 
above two oscillations, but the disk rotation also plays a part and gains or losses its energy. 
In the case of En ■ sign(u; — 7hQ)c > 0, the original oscillation is damped. Roughly speaking, 
in this case, both oscillations (original and intermediate ones) have the same signs in their 
wave energy. Hence, in this case, the main energy flow will be between disk rotation and the 
oscillations so that the oscillations are damped. It is noticed that in vertical resonance we 
always have En ■ sign(ci; — mf2)c > in the case where a resonant radius can appear in disk plane 
and the oscillations are always damped (see Kato 2004, 2008a). 

Finally, it is noted that the present resonant excitation mechanism has been applied to 
quasi-periodic oscillations observed in black-hole and neutron-star X-ray binaries, for example, 
by Kato and Fukue (2006), Kato (2008a) and Kato et al. (2008). 

The author thanks the referee for very careful reading of the manuscript with valuable 
comments, especially for pointing out widening of resonant region by pressure. 
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